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In doped lamellar cuprates, localized holes create ferro- 
magnetic bonds and cause spin canting similar to that caused 
by magnetic dipoles. At low temperatures, these dipoles 
are frozen, behaving like quenched correlated random fields. 
Renormalization group methods are used to show that such 
impurities cause a strong reduction of the two-dimensional 
antiferromagnetic correlations in the non-linear a model, and 
a related decrease in the three dimensional Neel temperature, 
in quantitative agreement with experiments. 
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Many lamellar copper oxides exhibit strong two dimen- 
sional (2D) Heisenberg antiferromagnetic (AFM) corre- 
lations, with a correlation length £^2D, and a 3D AFM 
ordering, at a Neel temperature T/y. Doping, which is 
believed to introduce a concentration x of electronic holes 
onto the oxygen ions in the Cu02 planes, leads to a rapid 
decrease in both ^2D and T/v, and to a disappearance of 
the AFM order above Xc ~ 2% This decrease has been at- 
tributed to frustration, due to strong ferromagnetic (FM) 
exchange on the Cu-O-Cu bonds which have localized 
holes . The FM bonds cause canting of the AFM mo- 
ments, which is described at long distance as resulting 
from effective dipoles |^J^ . This frustration was also pre- 
dicted to yield a magnetic spin glass phase for x > Xc, as 
recently confirmed in detail in doped La2Cu04 

At finite temperature T, the 2D AFM correlations of 
the undoped systems are excellently described by the 
classical non-linear-cr-model (NLcrM) with a renormal- 
ized stiffness constant ps |1C]. A few years ago, Glaz- 
man and loselevich (GI) [11| added dipolar defects to 
this model, integrated them out of the partition func- 
tion and described their effects on £^2D as a renormaliza- 
tion of ps- Although the locations of the dipole-like 
impurities and the unit vectors along the corresponding 
Cu-O-Cu bonds, a(r^), are randomly quenched, each im- 
purity involves an effective dipole moment Mm(r^) (with 
magnitude M) which is still free to reach annealed equi- 
librium in the presence of all the other dipoles. GI used 
an annealed averaging, expanded the renormalized ps to 
leading order in xps/T and used qualitative arguments 
to deduce a reentrant phase diagram for Tn{x). Here 
we argue that at sufficiently low T the dipole moments 
(which interact via randomly quenched dipole-dipole in- 
teractions) freeze in a random spin glassy way on the rel- 



evant length scales Therefore, they should also be 

treated as quenched variables. We thus solve the classi- 
cal NLcrM in the presence of such quenched random per- 
turbations, using the renormalization group (RG), and 
obtain S,2D for all the relevant values of T and x. In 
principle, one should use the annealed averaging for the 
first few RG iterations, until the renormalized cell con- 
tains more than a few impurities, and then (at low T) 
switch to the quenched averaging. To leading order in 
X, the two types of average give the same results, and 
our results agree with GI. However, at larger x our re- 
sults differ qualitatively from those of GI, and agree with 
experiments. 

Specifically, we express our results in terms of the 
renormalized parameters t = T/ps and A = Ax, where 
A is a dimensionless number. For t < A we find that 
S.2D{t, A) is independent of t, and is given by 



6D(i,A) »C(0,A)exp(27r/3A). 



(1) 



where C(0, A) « CoA'^. Fig. 1 shows this expression, 
with Co = 2.8A, lo — 0.8 and A = 20, in comparison 
with experimental and numerical results. 

At higher T, we find 



^2D{t,X)^C{t,X)exp(^^^ 
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with a slowly varying C(t, A). This agrees with GI to 
leading order in X/t, and its exponential part matches 
continuously with Eq. (|l|) at i = A. Fig. 2 shows this 
expression, using the known undoped values C {t, A) = 
1.92 A and 2t: ps = ISOmeV from Ref. |], for x = 
0.0225, 0.029 and 0.036, compared with experimental 
data for x = 0.020(5), 0.030(5) and 0.040(5) [|. The 
figure also shows our theoretical predictions for t < X, 
from Fig. 1 (the use of two different prefactors causes 
a discontinuity around t = X, smoothed by the dotted 
lines). The agreement with the data is comparable with 
(if not better than) that of the approximate expression 
(T, x) = C"^ (0, x) + (T, 0) used in Ref. [| with no 
theoretical basis (shown by the dashed lines). 

The 3D Neel temperature T]y{x) is roughly given by 
the solution of a^l^, = 1, where a represents either 
the relative interplane exchange, J±/J J±/2ttPs, or 
the in-plane spin relative anisotropy. Using a — 10~^ 
for La2Cu04 [§, Eq. (|l|) now yields x^ = 0.0183, in 
good agreement with the experimental value Xc = 0.0175 
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from Ref. |^ (but in conflict with Xc = 0.027 from Ref. 
1^). Approximating C{t,X) by the constant C(0, Ac) — 
(7(0,0.366) = 1.26A, the above expressions yield Tn{x) 
as shown in Fig. 3. Using the two different prefactors 
mentioned after Eqs. (|l|) and ^ would imply changes in 
T/v by less than 10%. Given that all the parameters were 
already determined from the experiments on ^2D, our 
T/v (a:) is in very good agreement with the measurements 
of Refs. 1^:^, based on Sr-doped La2Cu04. The disagree- 
ment with the data of Ref. may be due to possible 
systematic errors in x (these data were on 0-doped sam- 
ples, with X determined by Hall measurements). In con- 
trast with GI's heuristic arguments, which led to a reen- 
trance, our theory predicts that the line Tn{x) remains 
vertical at x = Xc up to the value where — Xc = Axc- 

The conventional method to treat the NLctM expands 
the order parameter unit vector n(r) of antiferromag- 
netism about a spatially uniform ordered state |p^ , |l7| . 
However, since the dipole-doped system has no AFM long 
range order in 2D even at T = 0, we preferred the orig- 
inal RG approach by Polyakov |l8|, which preserves ro- 
tational symmetry at all steps of the RG procedure. In 
fact, the conventional method generated local uncorre- 
lated random fields which should not be there [ p^ . 

Following GI, we describe the system by the reduced 
Hamiltonian 

= Wpurc + ^int, (3) 

where Tipurc is the classical NLcrM for the undoped sys- 
tem, representing the long wave length Hamiltonian re- 
lated to the fiuctuations of n(r), 

Wpurc-^ /'dr^(a,n^)2, (4) 

where i = 1, ...,c? and fi — \^...M run over the spatial 
cartesian components and over the spin components, re- 
spectively, and di = d/dxi. Hint was constructed by GI 
to reproduce the dipolar effects discussed in Ref. , 

nint^\ j dvY,Uv)-d,n, (5) 

with 

f,(r) = A/^(5(r-r,)a,(r,)m(r,). (6) 

We now decompose n(r) into a slowly varying part, 
given by the unit vector n(r), and (A/" — 1) fast variables 
(/)^(r) such that 

A/'-l 

n(r) = n(r)v/l-02(i.)+ ^ 0,,(r)e^(r), (7) 

with (jy^ir) ~ l]^j^0^(r). The unit vectors n(r) and 
e^(r), ji — 1, ...,7V — 1, form an orthonormal basis. The 



Fourier transform of the fast fields 4>^ is restricted to wave 
vectors q in the range < q < 1. The upper bound is 
the inverse of the microscopic length, and 6 = is the 
length rescale factor for the renormalization procedure. 
These q values are to be integrated out. After the itera- 
tion, the correlation length ^ is renormalized into 

The next stage involves integration over the fast vari- 
ables 4)^ [Q. This yields, hke in Ref. ||ll|, an effective 
2D dipole-dipole-type interaction between dipoles sepa- 
rated by a distance r < b. The interaction decays as r~'^. 
This renormalization can be repeated up to length scale 
b = Lq, and the resulting dipole-dipole interaction then 
applies to all the moments within the renormalized cell 
of size Lq. Clearly, to have a non-zero interaction within 
the cell we need Lq > x~^^'^, where x-i/'' is of the order 
of the average distance between impurities. We also re- 
quire Lq <^C S,2D, or else the impurities have no effect on 
the "pure" Heisenberg model behavior. 

A central issue is the averaging over the random fields 
f^. This includes averaging over the quenched random 
variables rg and a(r^), and over the variables m(r^), 
which are coupled via the above dipole-dipole inter- 
action. We use [ai(ri)aj(r2)] = Sij5{ri2)x/d, where 
[...] denotes the quenched average, ri2 = ri — r2, and 
find [ftf,{ri)fj^{r2)] = (5jj(5(ri2)A^^.(ri), with A^^(r) = 
mfj_{r)m,y{r)M'^x/d. At temperatures higher than the 
dipole-dipole interaction at mean distances (with en- 
ergy Edd of order Xps), one can follow GI and treat 
the m^'s as annealed. To leading order, this implies 
that (m^(ri)TO^(r2)) = S^i,S{ri2)/Af, and thus (A^^) = 
dfj,,yA = S^yM'^xj {dM). Given this relation, one derives 
the same leading order recursion relation for t as pre- 
sented below for the quenched case, while A remains un- 
renormalized in 2D [ p^ . 

We now argue that for sufficiently low temperatures 
one should treat the m's as quenched variables. It was 
argued (see Refs. and references therein) that spin 
glasses with long-range interactions, decaying as r""*, are 
at their lower critical dimension. Hence, the spin-glass 
correlation length, ^^g, increases exponentially with de- 
creasing T/Edd- Therefore, when T ^ Edd, we expect 
S,sg ^ Lq- When the cell size Lq is larger than x^^/'', 
each cell contains many impurities, and the strong spin- 
glass correlations will yield freezing of the dipole mo- 
ments (m(rfc)) within each cell into an apparent spin 
glass phase, with an Edwards- Anderson order parame- 
ter Q = [(m^(rfe))^]. At very low T, this becomes 
Q w 1/JV. We thus end up with a renormalized lattice 
in which the variables m(rfe) are frozen in random direc- 
tions, and all the factors in fj become quenched. Using 
the quenched average [m^(ri)m,y(r2)] — 5f^,j5{vi2)Q in 
A^jy, we now have 

[/v(i"i)/jv(r2)] = A(5^^(5ij(5(ri2), (8) 

where A = M^Qx/d « A. Fourier transforming the 
variables in Eq. (|5|), one can see that TLmt represents 
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random fields with quenched correlations of the form 
[h^ik)Kik')] = XkH^Jik + k'), where = 
Such correlations shift the lower critical dimension of 
the random field Heisenberg problem from 4 down to 
2 j2^j2l[], and at 2D one expects Imry-Ma domains of 
sizes given by Eq. (|l|) above. Indeed, our detailed RG 
calculations confirm this expectation. 

Integrating out the fast variables (f), and rescaling the 



momenta, yields the one-loop recursion relations |19 



dt 

dl 



A/"- 2 2 A/"- 1 , 
-et-\ r H tX, 



2n 



2n 



dl 2tt 27r 



(9) 



to leading order in e = d — 2. The rcnormalization proce- 
dure also generates many new random terms, which were 
not included in our initial H. However, all of these are 
irrelevant in the RG sense [ p^ . 

We now proceed to calculate ^2D- The prefacing an- 
nealed iterations up to length scale Lq = e^" are only 
expected to change t and A by factors of order unity. 
Thus, for £ > £o we solve Eqs. (^) with the approxi- 
mate initial values i(^o) = Iq = t -\- 0{t^,X) w t and 
A(^o) = Ao = Ax. The solution is particularly simple for 
the Heisenberg system, JV = 3. At 2D one finds 



A(^) = Ao( 1-^(^-4) 



m = m 



1 



Ao 

to 



1 



A(^)n^i/3 

~0 



(10) 



There is only one fixed point, at [X,t] = [0,0], and both 
quantities t{£) and A(^) flow away from it as £ increases. 
Using the standard scaling relation 

at,X)^e'mf^),X{i)), (11) 

we can find ^ from the matching condition 

max(t(r),A(r)) = 27r, (12) 

where £,{t{£*), X{£* j) = C{t,X) is a slowly varying func- 
tion of t and A. For our 2D solution we find 



/ 27r 



^2D{t,X) = LoC{t,x)exp{—{l - z^) 



(13) 



where — Ao/(27r) if Ao > to (i-e. X{£*) = 27r), and z is 
the solution of the equation 



+ (Ao/to - 1)2' - Ao/27 



(14) 



if Ao < tQ {t{£*) — 2tt). In the former case, this yields 
the exponential part in Eq. (^ (where we have dropped 
the subscript in Aq). From Eq. (p^, C{0,x) contains a 
factor Lq ^ x^^^'^. Another x-dependent factor will arise 
from higher order loops [|l^. Thus, C(0, x) ~ CqA'^, with 
uj of order unity (to be calculated or fitted). 



When Ao < to- so that z is not too small, one has z « 
1 — Ao/to, and thus our Eq. (|l3|) reduces to Eq. (||), where 
C{t,x) absorbs all the slowly varying additional factors. 
In this limit, the prefactor C{t,x) « C{t,0) is known 
from the two-loop and three- loop calculations. 
For La2Cu04, C(i,0) = 1.921. 

Finally, we discuss the derivation of the 3D Tn{x) in 
Fig. 3, as deduced from a^f^, ~ 1. The critical value Xc 
(above which there is no AFM long range order at any 
temperature) is given by the solution of 



A, = -47r/[31n(C'(0,A,)2a)], 



(15) 



with Ac — Axc- In fact, Eq. (|^) is expected to give a 
good approximation for ^20 for a range of temperatures 
obeying Ao > to. Therefore, the critical line t^ix) is 
expected to be practically vertical when t^^x) < Ax. 

For small x, i.e. Ao <S to, we can use Eq. (^, and 
obtain the linear approximation 



tN{x)/tN{0)^l-Ax/tN{0), 



(16) 



where tAr(O), the Neel temperature of the pure AFM, is 
the solution of tAr(O) = 47r/ ln[l/aC(tAr(0), 0)^]. 

For intermediate values of x (below Xc) Eq. (0) yields 



tNix) = X{l- 



i-^m 



aC{tisi{x), xy 



1/3- 



(17) 



The parameters used in Figs. 1-3 were chosen as fol- 
lows: The linear dependence of (t/27r) ln(^/C) on 1/t 
for a sample of La2Cu04+5, with = 90K gives 
A = 0.29 From Refs. ||,|l|l we deduce that this 

value of Tat corresponds to a; w 0.0145, hence A « 20. 
This is the only free parameter used to describe both ^20 
and Tn{x) for t > X\ The values of Co and ui were chosen 
to fit the data in Fig. 1. 

At temperatures lower than ^ 200 — 250K, £^20 does 
not depend on the temperature up to exponentially small 
terms, of the order of ^(0, T) ~ ^ [^j . This nontrivial prop- 
erty is reproduced well by the theory. At the highest tem- 
peratures the calculated values of £20 are smaller than 
the measured ones. Perhaps, at such high temperatures 
thermal fluctuations, which should decrease the dipole 
moment, become important. In fact, at these tempera- 
tures one might have to allow some annealed averaging, 
both on the dipole moments and on their locations. 

In conclusion, we employed the NLctM to treat the 
effect of hole doping on the properties of lamellar copper 
oxides. Our theory seems to agree with experiments on 
both £2D{t,X) and Tn{x). 

We have benefitted from many discussions with 
R. J. Birgeneau, A. B. Harris, A. S. loselevich, 
M. A. Kastner, D. E. Khmelnitzkii, V. V. Lebedev, 
T. Nattermann and M. Schwartz. This project has been 
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FIG. 3. Tn{x)/Tn(0) versus x. Full line is theory, and 
the points are from experiments: full circles from Ref. 4 , -|-'s 
from Ref. 5, empty circles from Ref. 15. Dashed line indicates 

Ao = to- 
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FIG. 1. Dependence of (,2d at T = on x. The empty[2] 
and full [13] circles indicate experiments, +'s show numerical 
simulation[14] data. The solid line represents Eq. (1) with 
C(0, A) = 2.8A° *A and A = 20a;. 



FIG. 2. Depeudeuce of ^2£)(t, A) on T for several concen- 
trations. Symbols are from experiments: empty circles for 
X = 0.02, full circles for x = 0.03, +'s for x = 0.04, all from 
Ref. 2. Full lines show Eqs. (1) (with C = 2.8A'' *A) and (2) 
(with C = 1.92 A), and dotted lines interpolate between these 
low and high temperature theories. 
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